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We propose a hybrid model for medium-induced parton energy loss, in which the hard scales in 
the process are treated perturbatively, while the soft scales which involve strong coupling dynamics 
are modeled by AdS/CFT calculations. After fitting a single parameter on Raa for central Au+Au 
collisions, we are able to predict different observables like Raa and Iaa as a function of centrality and 
reaction plane. We obtain a consistent picture of how jet quenching is modified if the quark-gluon 
plasma is strongly interacting, and we provide quantitative predictions. 



I. INTRODUCTION 

o \ 

Understanding the quark-gluon plasma (QGP) created in nucleus-nucleus collisions [lHl] at the Relativistic Heavy 
Ion Collider (RHIC) is a challenging task. Many measurements have been carried out to provide insight on the 
properties of that dense QCD matter, yet a few years after the discovery, the question whether the QGP formed at 
P-h • RHIC is weakly or strongly coupled remains unanswered. On the one hand, some features of the plasma are naturally 
' understood from strong cross-sections, and even fit quantitatively within a perfect-liquid picture, like the large elliptic 
flow On the other hand, one would like to establish a global picture and be able to also answer the question 

using hard probes. While qualitatively the medium is opaque to energetic particles, implying strong interactions, a 
quantitative strong-coupling description of jet quenching is missing. 

Hard probes are believed to be ideal processes to study the properties of the QGP, they are thought to be understood 
Oh| well enough to provide clean measurements. Observables built to measure medium effects on particle production, like 
the nuclear modification factors Raa, are not always easily reproduced in perturbative QCD (pQCD) calculations, 
even for hard production. For instance in older computations, to reproduce light-hadron Raa's within the pQCD 
framework of medium-induced energy loss the value of the so-called jet quenching parameter q is adjusted to 

■ 5—10 GeV 2 /fm [llj. It is otherwise estimated to be smaller (1 — 3 GeV 2 /fm) for a weakly-coupled plasma at RHIC 
temperatures. A more recent comparison study of different pQCD approaches within the same medium model fl2j ] 
shows that assumptions about how q depends on thermodynamical parameters play a crucial role and that in principle 
pQCD based models may also work with a low q. While our current understanding of the pQCD picture must be still 
be improved , it is unclear if the pQCD approach can describe the suppression of high— Pt particles at RHIC. 

The goal of this letter is to quantify the expected modifications of the fully-perturbative radiative energy loss 
results in QCD, if the plasma is strongly coupled. Addressing the strong coupling dynamics in QCD is an outstanding 
problem, lattice simulations remain the only method to obtain quantitative results, but are inefficient when real- 
time dynamics is required, which is the case when an alyzing jet quenching. However for a class of non-abelian 
thermal gauge theories, the AdS/CFT correspondence |13T - [l5| provides an alternative: it allows to investigate the 
strong coupling regime in the large- N c limit, essentially by mapping the quantum dynamics of the gauge theory into 
analytically-tractable classical gravity dynamics in the fifth dimension of a curved anti-deSitter (AdS) spacetime. 

We propose a hybrid model for medium-induced parton energy loss, in which the hard scales in the process are 
treated perturbatively, as in the standard pQCD radiative energy loss calculations, while the interaction with the 
plasma which involves strong-coupling dynamics is modeled by AdS/CFT calculations, for the Af = 4 super- Yang- 
Mills (SYM) theory. While this theory is quite different from QCD (it is highly supersymmetric and conformal), it 
has become a popular approach to assume that the SYM plasma approximates well the QCD plasma, just above 
the critical temperature T c , where the conformal anomaly in QCD is small. Information about real-time dynamics 
within a thermal background can be obtained with this setup. The best-known example is the calculation of the shear 
viscosity to entropy ratio [lril [l7| , which in the strong coupling limit holds for any gauge theory with a gravity dual. 
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After a bare parton is created in a hard process, it starts to build its wave function with quantum fluctuations. 
The presence of the medium prevents the parton to become fully dressed as it would in the vacuum: long-range 
fluctuations, typically with virtuality less than the plasma temperature T, are screened out of the wave function. The 
parent parton then loses the energy that those fluctuations carry into the medium. We would like to compute the 
properties of this medium effect in QCD. In the high parton energy limit, the energy is transfered into the medium from 
the perturbative part of the parent parton wave function. Then we explicitly know what the quantum fluctuations 
are, quarks and gluons, and their dynamics can be computed from first principles 18]. On the contrary, if one were to 
probe the soft part of the wave function, then we wouldn't know how to deal with the fluctuations and their dynamics 
in QCD. 

Having in mind the description of high— Pt hadrons in heavy-ion collisions, which come from the fragmentation of 
energetic partons, the authors in [l9| proposed to use the Baier-Dokshitzer-Mueller-Peigne-Schiff-Zakharov (BDMPS- 
Z) energy loss formalism [2(| [2l| to calculate the gluon emission process, which involves a s {p T ), while describing the 
interaction of the emitted gluons with the strongly-coupled medium with AdS/CFT calculations. Their motivation is 
that, in the high-parton energy limit, the BDMPS-Z analysis of medium induced parton energy loss is under control, 
and whether the plasma is weakly-coupled or not gluon radiation is the dominant mechanism |22|. Moreover, applying 
strong-coupling techniques to the entire radiation process would be inappropriate, because QCD is asymptotically free, 
and the calculation of the gluon distribution in the parent parton wave function is correctly treated with perturbative 
methods. 

In the so-called multiple soft scattering approximation of the BDMPS-Z calculation, the properties of the energy loss 
are expressed in terms of the transport coefficient q, a constant rate at which gluons in the parent parton wave function 
pick up transverse momentum squared dp 2 ± /dt — q cx T 3 . Instead of evaluating q perturbatively, the authors of [l9| 
proposed to evaluate this coefficient at str ong coupling. Unfortunately, at strong coupling, the energy loss cannot be 
expressed in terms of a constant dp 2 j_/dt [23l |2~3|. However, an interpretation based on the multiple soft scattering 
approximation in the BDMPS-Z calculation is still possible, but instead with multiple scatterings characterized by the 
constant rate d\p±\/dt oc T 2 . Our strategy is to express the properties of the energy loss in terms of the accumulated 
p\ as in the BDMPS-Z calculation, and then to evaluate this accumulated transverse momentum squared, that we 
shall denote Q 2 , at strong coupling: Q 2 = T 4 L 2 (instead of Q 2 = qL at weak coupling), where L is the extent of the 
medium. Our hybrid model is then completed by a realistic account of the plasma geometry and expansion, as in 25] . 

The plan of the letter is as follows. In Section II, we recall the different pieces of the BDMPS-Z medium-induced 
radiative energy loss calculation, as well as the physical picture emerging from such a weak-coupling analysis. We 
also explain what part of the calculation should be modified in the case of a strongly-coupled plasma, and propose 
a substitution inspired by AdS/CFT results. In Section III, we discuss how the plasma geometry and expansion are 
implemented in order to obtain meaningful predictions. In Section IV, we present results for Raa and Iaa at high 
Pt as a function of the reaction plane, observables for which the stronger path-length dependence of the energy loss 
at strong coupling, compared to weak coupling, produces sizable differences. Section V is devoted to conclusions and 
outlook. 



II. MEDIUM-INDUCED PARTON ENERGY LOSS 



We denote the energy of the hard parton E, and work in the large parton energy limit. We also denote oj and 
k± , the energy and transverse momentum of the virtual gluons in the parent parton wave function (note that we are 
using pt to denote transverse momenta with respect to the collision axis, and p± with respect to the jet axis). What 
prevents the radiated gluons to have an arbitrarily large energy, is that their coherence time t c should be smaller than 
the length of medium L the parent parton goes through before exiting the medium: 

t c = uj/k'i , uj < L k]_ . (1) 

The transverse momentum squared, or virtuality, of the radiated gluons is also limited: the interaction with the plasma 
cannot overcome an arbitrary large virtuality and put on shell any fluctuations. Let us denote p^_(i c ) the transverse 
momentum squared acquired by a gluon of coherence time t c , then only the gluons with virtuality k\ < p^_(i c ) can 
be radiated. For those emitted gluons which dominate the energy loss of the parent parton, we shall denote their 
virtuality and energy 

Q 2 s =pl(L), and u m =LQ 2 s . (2) 

The saturation momentum Q s depends on the temperature T and length L of the plasma, and the precise way it 
depends on it is different if the plasma is weakly coupled or strongly coupled. 
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The first part of the BDMPS calculation is to determine the energy distribution of the radiated gluons cudl/dus, 
from which the total energy loss AE, as well as the quenching weight 

P(AE) = £ - 
n! 

71=0 

necessary to compute Raa, can be computed [26|. The saturation momentum is what determines the properties of the 
energy loss, as the quenching weight are expressed in terms of Q 2 and oj m . The other part of the BDMPS calculation 
assumes that the plasma is made of thermalized weakly-interacting quarks and gluons, which allows to determine Q 2 
and uj m as functions of the plasma temperature and length. 




A. In a weakly-coupled plasma 

At weak-coupling, medium- induced gluon radiation is due to multiple scatterings of the virtual gluons [20] : if, while 
undergoing multiple scattering, the virtual gluons pick up enough transverse momentum to be put on shell, they 
become emitted radiation. The accumulated transverse momentum squared picked up by a gluon of coherence time 
t c is 

pi(t c ) = fi 2 j = q t c , which gives Q 2 S = qL , (4) 

where fi 2 is the transverse momentum squared picked up in each scattering, and A is the mean free path. These 
medium properties are involved through the ratio q = /i 2 /A, this is the only relevant information about the medium. 
In terms of the temperature T, one has q ~ ct s T 3 . The total energy lost by the quark is therefore AE oc a s qL 2 , where 
the factor a s comes from the probability for the gluonic fluctuation in the wave function, and essentially counts the 
number of gluons which lose the energy qL 2 . 

An important quantity which enters the calculation of udl/dio is the scattering amplitude of a color singlet dipole 
off the plasma, which at weak coupling is given by 

N(r) = 1 - e^ Lr2 . (5) 

This explains why we denoted the accumulated p\ by Q 2 : this sets the scale at which the scattering becomes strong. 
The dipole is made of the parent parton and its antiparticle (it could be a qq dipole or a gg dipole, the difference is 
a color factor absorbed in q). In the so-called multiple soft scattering approximation that we are considering here, 
the scattering amplitude for small dipoles in a single scattering N(r) ~ qXr 2 — fi 2 r 2 is explicitly used. This allows 
to obtain the quenching weight as a function of Q 2 and u> m [27] . In phenomenological calculations, the geometry and 
expansion of the plasma are taken into account by integrating over the parton trajectory £ in the medium [251 ]: 

poo p oo 

Q 2 (v , 4>) = K I dt T 3 (0 , w m (r , 4>) = K / d£ £ T 3 (0 . (6) 
Jo Jo 

Tq and <p denote the position of the hard process that created the parton and its direction with respect to the reaction 
plane. The normalization K is the only parameter in the problem. 



B. In a strongly-coupled plasma 

We would like to address the case of a strongly-interacting quark-gluon plasma. We propose to model the strong 
coupling dynamics by replacing in the BDMPS calculation, the values of Q 2 and uj m by their strong-coupling expres- 
sions. In other words, we assume that what is radiated into the medium still comes from the perturbative part of the 
energetic parton wave function, but we assume that the interaction of those virtual gluons with the medium, and the 
way they are freed, is governed by strong-coupling dynamics. At strong coupling the value of Q s can be estimated 
by a classical gravity calculation, in 5 dimensions. Indeed, using the AdS/CFT correspondence (l3l - [T5j | which maps 
the quantum dynamics of the SYM field theory onto classical dynamics in a fifth dimension, one can evaluate the 
saturation momentum for this theory. We shall then use it as an input for our model. In those classical calculations, 
the plasma acts as a force of magnitude T 2 on the energetic parton and its wave function, to which correspond a 
string in the fifth dimension. 
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The work done by the plasma on a fluctuation of coherence time t c is T 2 t c . To put a quantum fluctuation on shell, 
this work should overcome its virtuality, therefore we identify 

pi (t c ) = TH\ , which gives Q\ = T 4 L 2 (7) 

for those fluctuations which dominate the energy loss. The expression p\ (t c ) was obtained in [28[ with an R— current 
deep inelastic scattering (DIS) calculation. How this result translated into the length dependence of the saturation 
scale and the energy loss was discussed in [23|,[24| for an energetic heavy quark, and the main result was that at strong 
coupling AE oc L 3 instead of the L 2 law at weak coupling. Calculations of DIS off a finite-length strongly-coupled 
SYM plasma confirmed the result (29W31I ] . It is obtained that Q 2 = T 3 L/x, where x is the Bjorken variable. In DIS 
off a weakly-coupled QCD plasma, the result Q 2 = T 3 L is independent of x, i.e. independent of the collision energy 
\fs ~ 1/x. This comes from the fact that at weak-coupling, the quasiparticle which determines high-energy scattering 
is a spin 1 object (usually called Pomeron), and therefore s J_1 ~ 1 (there will be an x dependence via the plasma gluon 
distribution xg{x) if one takes into account quantum evolution |33|1. By contrast, at strong coupling the quasiparticle 
which determines high-energy scattering is a spin 2 object (therefore called graviton), and s" 7-1 ~ which explains 
the differences between the weak- and strong-coupling saturation momenta. In the energy loss calculation, the extra 
1 jx factor is simply replaced by t c T, and for a plasma of small enough extent such that it is indeed the finite length 
of the plasma which limits the coherence time of the freed fluctuations, one recovers Q s = T 2 L. Finally, the AE oc L 3 
law at strong coupling can also be obtained for gluons and light quarks, from the dynamics of falling strings (33l. |34|. 

In order to make the connection with results in the literature, let us briefly consider the case of an infinite-extent 
plasma. Then t c = uijk\ is not restricted anymore, and the transverse momentum bounds become k\_ < qui in the 
weak coupling case, and kj_ < T 2 uj in the strong coupling case. These, coupled with the wave function suppression 
for uj > Ek±/M (the dead cone effect) [35| in the case of a heavy quark parent parton give Q 2 = (qE /M) 2 / 3 at weak 
coupling and Q 2 = T 2 E/M at strong coupling for the virtualities of the most energetic freed fluctuations, whose 
energies are ui m — EQ S /M. If the plasma has a finite length, but bigger than the corresponding coherence times 
Ej (MQ S ), then the matter is effectively of infinite extent. In the strong coupling case, we recover that the saturation 
scale Q s — T^jE/M in the gauge theory corresponds on the gravity side to the st ring world sheet horizon of the 
trailing string calculation [36l l37| , which determines the heavy quark energy loss [38], l39{ . 

When computing dl/dui, we are still using the multiple soft scattering approximation, but with the strong coupling 
expression for Q 2 . At the end, we substitute into the analysis of [25l |. the following expressions: 

/•OO /"OO 

Q 2 (v , 4>) = K / d£ £ T 4 (£) , Lj m (r , cf>) =K d^ 2 T 4 (£) . (8) 
Jo Jo 

As done in standard weak-coupling calculations, the normalization will be fitted on Raa data as a function of 
transverse hadron momentum Pt for 10% central 200 AGeV Au-Au collisions for one value of Pt- This value of K 
is then used to predict the Pt dependence, with a different hydrodynamical run the centrality dependence and also 
other observables such as back-to-back correlations. 

Let us make a final comment. At weak coupling, gluons accumulate transverse momentum via a multiple scattering 
process characterized by the constant rate dp^/dt = q ~ T 3 , independent of the length of the medium (they receive 
random p± kicks). At strong coupling, the (faster) accumulation of transverse momentum can also be viewed as 
a multiple scattering process, but characterized by the constant rate d\p±\/dt — T 2 (gluons receive |pj_| kicks). 
Therefore, dp\/dt ~ T 2 |p^| (— > T 2 Q S for the dominant fluctuations) is not a constant and should not be thought of 
a strong coupling result for the jet quenching parameter q. For a finite- length plasma one has dp\ /dt ex. T A L while 
for an infinite- length plasma the result is dp\jdt oc T 3 ^J 1 EjM (we recover the results of (4p| - l42l | ) . Even though for 
infinite-extent matter the result looks similar to the weak coupling result, the physics is different. 

III. IMPLEMENTING THE PLASMA GEOMETRY AND EXPANSION 

In order to compute observables from the quenching weight Eq. @, we have to evaluate this expression for any 
given path of a parton through the medium using the line integrals in Eqs. ([5]). average over the paths through the 
medium characteristic for a particular high— Pt observable and finally fold the averaged energy loss probability with 
the QCD expressions for hard hadron production. In the following, the medium evolution is always specified in terms 
of a 3-d hydrodynamical model [43| which is well constrained by a large number of bulk observables. From this model, 
the local medium temperature T at any spacetime point can be determined for Eqs. (|8"|). For comparisons, we will 
also use Eqs. ©• In this case T will be substituted by e 1 / 4 where e is the energy density to make the exact connection 
to previously published results where this scaling was assumed. This is known to make a difference in terms of fitting 



K but it does not change the shape of the resulting suppression curves, see [12[ for a detailed discussion 
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The computation is for technical reasons different for single hadron suppression in terms of the nuclear modification 
factor Raa and for dihadron correlations in terms of the suppression factor Iaa ■ While the former quantity is obtained 
with standard multi-dimensional numerical integration techniques, the latter quantity is solved in a Monte Carlo (MC) 
simulation. A detailed description of the model is found in [25|], here we outline the essential steps. 

From the quenching weight P(AE) given a single path in the medium as computed using Eqs. (J5J, we define the 
averaged energy loss probability distribution (P(AE))t aa as 

1 />27T />OG />00 

(P(AE)) Taa =— / d<\> \ dx / dy P(x Q ,yo)P(AE). (9) 

Z7T JO J-oo J-oo 

The averaging weight P{x$, yo), i.e. the probability for finding a hard vertex at the transverse position ro = (xo,ya) 
and impact parameter b, is given by the product of the nuclear profile functions: 

p( , T A (v +h/2)T A (r -h/2) 

P(X0 ' y ° ) = IWbl ' (10) 

where the thickness function is given in terms of Woods-Saxon the nuclear density pa(t, z ) as T A (r) — J dzpA(r, z). 
The energy loss probability P(AE) is derived in the limit of infinite parton energy. In order to account for the finite 
energy of the partons we truncate (P(AE))t aa at AE = Ej Ct and add <5(A — -Ejot) deP(e) to the truncated 
distribution to ensure proper normalization. The physics meaning of this correction is that we consider all partons as 
absorbed whose energy loss is formally larger than their initial energy. We calculate the momentum spectrum of hard 
partons in leading order perturbative QCD (LO pQCD) (explicit expressions are given in [25| and references therein). 
The medium-modified perturbative production of hadrons can then be computed from the expression 

da^ h+x =J2^lf^ f+X ® (P(AE)) Taa ® Dyc h (z,n%) (11) 
/ 

with Dj°^ h (z, /ip) the fragmentation function with momentum fraction z at scale p? F [44j], and from this we compute 
the nuclear modification factor Raa which is the ratio of hard hadron production in A-A collisions normalized to the 
production in p-p collisions scaled by the number of binary collisions, i.e. 

dN h AA /dP T dy 

Raa{Pt ^ )= T A A(b)da PP /dP T dy (12) 

For the computation of in-medium back-to-back correlations with a hard triggered hadron, the trigger bias has to be 
included consistently into the averaging over geometry. For this purpose, we simulate a large number of back-to-back 
events in the medium. We start with the pQCD expression for the production of partons k, I from two colliding 
objects A, B 

, 2 , J = 2^X 1 f i/ A{x 1 ,Q 2 )x 2 f 3 /B(x2,Q 2 ) -J (13) 

dp^dyidy 2 ^ dt 

where fi/A{ x i,Q 2 ) is the parton distribution function at fractional momentum x\ and scale Q 2 which is different in 

the proton [i^, |46| and in the nucleus [U |48| and — are perturbatively calculable cross-sections for the various 
hard partonic subprocesses. We MC sample this expression at midrapidity y\ = y 2 = to get a back-to-back parton 
pair with given parton type and momentum px- At the same time, the spatial position of the event is MC sampled 
from Eq. pH|) . We propagate each of the partons through the medium and compute P(AE) using Eqs.©, then 
determine probabilistically a value of AE from P(AE) and subtract it from the parton energy before hadronization. 
If AE > E, i.e. if the parton energy loss is larger than its energy, we consider the parton absorbed by the medium, 
otherwise it emerges with a finite energy and we hadronize it subsequently. 

In order to determine if there is a trigger hadron above a given threshold, given a parton k with momentum px, 
we need to sample A\~^ h {zi 1 pT), i.e. the probability distribution to find a hadron h from the parton k where h is the 
most energetic hadron of the shower and carries the momentum Pt = Zi • Vt- We extract A\(zi,pt) from the shower 
evolution code HERWIG [4!|. The procedure is described in detail in [5(3]. Sampling A\(z\,pr) for any parton which 
emerged with sufficient energy from the medium provides the energy of the two most energetic hadrons on both sides 
of the event. The harder of these two defines the near side. The hadron opposite to it is then the leading away side 
hadron. In order to compute the correlation strength associated with subleading fragmentation of a parton emerging 
from the medium we evaluate A 2 (zi, z 2 ,Pt) (also extracted from HERWIG), the conditional probability to find the 
second most energetic hadron at momentum fraction z 2 given that the most energetic hadron was found with fraction 
z\. Given a trigger, we can then count the number of events in which correlated hadrons on near and away side fall 
into certain momentum windows and determine the per-trigger yield. The dihadron suppression factor Iaa is then 
computed by dividing the per-trigger yield in A-A collisions by the per-trigger yield in p-p collisions. 
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• PHENIX data 

— ASW/AdS, b = 2.4 fm, e 3/4 /T 4 scaling 

— ASW, b = 7.5 fm, e 3 scaling, in plane 
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FIG. 1: Left panel: Nuclear modification factor Raa for 200 AGeV Au+Au collisions as a function of Pt, for central collisions 
(data and b = 2.4 fm) and 25% peripheral collisions (b = 7.5 fm). Results for the weak coupling (labeled ASW) and strong 
coupling (labeled AdS) scenarios for the medium are displayed (see text). The data are from [5lfl . Right panel: the same strong 
coupling predictions are compared with data for peripheral collisions and two reaction plane angles. The data are from [521. |53|. 



IV. RESULTS 



We show the Pt dependence of the nuclear modification factor Raa for 200 AGeV Au+Au collisions as compared 
with the PHENIX data for 10% central collisions [5l[ in Fig. [1] (left plot) for both the weak coupling and the strong 
coupling assumption for the medium. In both cases, with a single fit parameter, a good (and virtually identical) 
description of the data can be made. However, in the extrapolation to more peripheral collisions (impact parameter 
b = 7.5 fm), differences between the two scenarios become apparent, dependent on the angular orientation of the 
high— Pt hadron with respect to the reaction plane. In particular, the spread between in-plane and out-of-plane 
emissions is increased in the strong coupling case, which is in agreement with preliminary data |52| . |53| as is shown in 
Fig. ED (right plot). 

This is consistent with the different pathlength dependence of energy loss which is for a constant medium L 2 in 
the weak coupling case but L 3 in the strong coupling case. Consequently, the average pathlength difference between 
in-plane and out-of-plane emission is weighted with a different power also in the case of the expanding hydrodynamical 
medium, which translates into a stronger effect in the strong coupling case. Because both calculations are adjusted on 
central collisions, they predict a similar suppression for out-of-plane emissions (the Raa's are similar to each other but 
both bigger than for central collisions where the temperature is higher), as these involve the same average pathlength. 
Consequently, the stronger L dependence causes the strong coupling scenario to predict less suppression than the 
weak-coupling scenario for in-plane emissions, where the pathlength is on average shorter. 

Based on the results of (55|, back-to-back correlations where the difference in average pathlength is maximized due 
to the surface bias of the trigger hadron would be an appropriate observable to show the differences between the weak 
and strong coupling scenario for the medium more clearly. The away side per-trigger yield for an 8+ GeV hadron 
trigger as measured in 200 AGeV central Au-Au collisions by the STAR collaboration [56| is shown in the left panel of 
Fig. [2] and compared with both the weak coupling and the strong coupling treatment of the medium. As expected, the 
different weighting of the pathlength difference in the two approaches translates into a stronger suppression for the 
L 3 pathlength dependence of the strong coupling approach, the net effect however is mild, especially in the region of 
associate hadrons at 6+ GeV momentum where fragmentation can safely be assumed to dominate hadron production, 
and also where our approach is best applicable. 

In the right panel of Figf2J we present a more differential observable: the away side dihadron suppression factor 
Iaa for various collision centralities and away side momentum bins, as calculated in (57j within the weak coupling 
approach to medium properties, and compared with the strong coupling approach of this letter. The presentation 
within a bin is arranged in two groups, the weak coupling results first, the strong-coupling results next. In each 
group, the color coding refers to the centrality of the hydro run (there is more suppression for more central collisions), 
and for each centrality there are two points (with different symbol shape) - the in-plane result is always the upper 
point, showing less suppression, the out of plane result is the lower point. The in-plane and out-of-plane results are 
closer to each other for more central collisions. Finally, one observes that the strong coupling approach predicts more 
suppression than the weak coupling one for central collisions (consistent with the per-trigger yield result presented 
above) but less suppression for peripheral collisions, and a larger spread between in-plane and out-of-plane emissions. 
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FIG. 2: Left panel: Per trigger yield on the away side with respect to a 8+ GeV hadron trigger in back-to-back correlations 
in central 200 AGeV Au+Au collisions for two different momentum bins. Shown are STAR data the calculated result for 
a weak coupling and for a strong coupling treatment of the medium. Right panel: Nuclear dihadron away side suppression 
factor Iaa for a 12-20 GeV trigger in 200 AGeV Au+Au collisions for different impact parameter and in-plane vs. out-of-plane 
emission as a function of away-side momentum bins. Shown are weak coupling (ASW) and strong coupling (AdS) treatment 
of the medium. The points have been shifted along the rr-direction within each momentum bin for clarity of presentation. 



For completeness, the values of the parameter adjusted to describe Raa for central collisions are K ~ 27 in the 
weak-coupling case (Eqs. ([6])) and K ~ 6 in the strong-coupling case (Eqs. ([8])). The weak-coupling value is a factor 
3.5 bigger than expected for an ideal gas [54| ; we are not aware of a formula with which we could compare the 
strong-coupling result, but the fact that K is of order one is satisfactory. 



V. CONCLUSIONS 



In this letter, we propose a new description of the quenching of high— px particles moving in a strongly-interacting 
QGP. In our model, energetic partons lose energy through the emission of long-lived gluons (with coherence time 
t c ~ L, the plasma extent), as dictated by perturbative calculations, while the amount of energy that these gluons 
carry away in the medium is enhanced by the strongly-coupled nature of the plasma, compared to weak-coupling 
calculations. We do not perform a direct AdS/CFT calculation, we rather follow the idea of Ref. [ljj, that whether 
the plasma is weakly-coupled or not, gluon radiation is the dominant mechanism by which high— px particles lose 
energy. However, the assumption in [19( to work with a constant q computed at strong-coupling does not allow to 
implement the path-length dependence of the energy loss, which is important for phenomenology. Our prescription 
does. 

With respect to a weakly-coupled plasma, the path length dependence of the energy loss of energetic partons in a 
strongly-coupled plasma is stronger (oj m — T 4 L 3 vs. ui m — qL 2 ). In order to quantify the effect that this has on physical 
observables like Raa and Iaa at high Pt, we have proposed to substitute Eqs.© in a successful phenomenological 
analysis at weak coupling, by their strong-coupling expressions Eqs.©. In both cases we have assumed a s (p T ) <C 1, 
the high-parton energy limit, and used the BDMPS-Z radiative energy loss formalism. Only for low— px observables 
can one hope to apply strong-couplin g te chniques to the entire process, and interesting results for dihadron correlations 
have also been obtained in this case [58l |59| . 

Our two calculations are fixed by one-parameter fits to Raa for central Au+Au collisions, resulting in good and 
comparable descriptions of the data. Differences between the two cases are then seen when considering non-central 
collisions: the spread between the in-plane and out-of-plane Raa's is increased in the strong coupling case, which is 
in good agreement with the data. The suppression of Iaa is also increased for central collisions, as well as the spread 
between the in-plane and out-of-plane values for non-central collisions. Our hybrid model allows to make quantitative 
predictions which could be tested with future measurements. One concern is that the theoretical uncertainties could 
be large enough to blur the potential resolving power of the proposed measurements, for which the quantitative 
differences between the two approaches are small. Moreover, since a major source of uncertainty in jet quenching 
calculations (and more globally in heavy- ion collisions) is the shape of the initial hydro profile (CGC vs. Glauber), 
quantifying the theoretical uncertainty is not straightforward. 
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In the future, it would also be interesting to investigate how our substitution impacts the Raa of heavy hadrons 
(or rather non-photonic electrons), which is difficult to accommodate in the pQCD framework, and one could also 
compare with results obtained fully from strong-coupling dynamics 60] . 
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